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ISOTHERMAL ELASTOHYDRODYNAMIC LUBRICATION OF POINT CONTACTS 
II - ELLIPTICITY PARAMETER RESULTS 


by Bernard J. Hamrock and Duncan Dowson* 

Lewis Research Center 

SUMMARY 

A numerical solution of the isothermal elastohydrodynamic problem for point con- 
tacts is presented which reproduces all the essential features of the previously reported 
experimental observations based upon optical interferometry. In particular, the two 
side lobes, in which minimum film thickness regions occur, emerge in the theoretical 
solutions. The influence of the ellipticity parameter on solutions to the point contact 
problem is explored. The ellipticity parameter k was varied from 1 (a ball on a plate) 
to 8 (a configuration approaching line contact). It is shown that the minimum film thick- 
ness can be related to the well known line contact solutions by a remarkably simple ex- 
pression involving either k or the effective radius of curvature ratio Ry/R^. 


INTRODUCTION 

Most of the work to date in elastohydrodynamic lubrication (EHL) has dealt with 
line contacts. Grubin (ref. 1) was the first to attempt a solution to the isothermal EHL 
line contact problem. In his analysis it was assumed that the shape of the elastically 
deformed surfaces was the same as the shape produced in a dry (Hertzian) contact. 

This assumption facilitated the solution of the Reynolds' equation in the inlet region to 
the contact and enabled the film thickness in the central region to be determined with 
acceptable accuracy. In references 2 and 3 an empirical formula for the solution to the 
isothermal line contact EHL problem was obtained. This formula shows the effect of 
speed, load, and material properties on minimum film thickness and is based on their 
theoretical solutions. The experimental observations concur with the minimum film 
thickness formula. 

The study of point contacts has mostly foUowed experimental lines (e. g. , refs. 


♦Professor of Mechanical Engineering, Leeds University, Leeds, Er^land. 


4 to 6). Before Archard and Cowking (ref. 7) presented their central film thickness 
formulation for point contact, it was common practice to use line contact minimiun film 
thickness to determine the film thickness in point contact situations such as in ball bear- 
ings and gears. Archard and Cowking adopted an approach similar to that used by 
Grubin (ref. 1) for line contact; that is, a Hertzian deformation is assumed when the 
shape of the surfaces is identical to that which occurs under the same load in the ab- 
sence of a lubricant. The film shape resulting from this t3^e of analysis is such that 
the central region is assumed to be parallel to the contacting plane, and the analysis is 
primarily concerned with deriving an expression for the central film thickness. Cheng 
(ref. 8) also used a Grubin approach in his point contact analysis. 

In an earlier report (ref. 9) we outlined the complete elastohydrodynamic lubrica- 
tion (EHL) point contact problem. This called for the simultaneous solution of the elas- 
ticity and Reynolds' equations. In the elasticity analysis the contact zone is divided into 
equal rectangular areas. It was assumed that a uniform pressure was applied over each 
element. In the numerical analysis of the Reynolds' equation, a ^ substitution, where 
(p is equal to the pressure times the film thickness to the three halves power, was used 
to help the relaxation process. Only the theory was presented, with the results to fol- 
low in later reports. 

The present report presents the first of the results. Herein the ellipticity param- 
eter is varied from 1 (a ball and plate curvature combination) to 8 (approaching line con- 
tact) while keeping the other dimensionless parameters fixed. Contour plots of pressure 
and film thickness are presented. From the results the minimum film thickness formula 
will be given relating the effect of the elliplicity parameter on film thickness. By using 
these results, a designer will have at his disposal a complete contour map of film thick- 
ness. Comparisons will also be made with Archard and Cowking (ref. 7). 


SYMBOLS 


a 

b 

E 

E' 


F 

G 

H 


semimajor axis of contact ellipse 
semiminor axis of contact ellipse 
modulus of elasticity 


V{[(‘ ■ "i)K] +[(i - 



normal applied load 

dimensionless material parameter, E'/p.^ 
dimensionless film thickness, h/R 


2 


H; 


H, 


TT 

min 

TJ . 

min, L 
h 

k 

P 


P . P . 
mm, min 


P 

P 

R 

r 

U 

u 


iv, as 


dimensionless central film thickness obtained from ref, 7 

dimensionless central film thickness for point contact 

dimensionless minimmn film thickness for point contact 

dimensionless minimum film thickness for line contact 

film thickness 

ellipticity parameter, a/b 

dimensionless pressure, p/E' 

defined by eqs. (20) and (22) 

pressure 

asymptotic isoviscous pressure 
effective radius 
radius of curvature 

dimensionless speed parameter, Vp^RE' 
surface velocity in x direction 


V 

V 

W 

x, X, X,1 

y, Y, Y J 


+ v^ 

surface velocity in the y direction 
dimensionless load parameter, F/E’Rj^Ry 

coordinate systems used in the report 


Z 

01 ,^ 

9 

A 

V 

Pq 

Subscripts: 


viscosity pressure index, a dimensionless constant 
constants used to define the density 
atmospheric viscosity 
spin to roll ratio defined in eq. (8) 

1/[1 + (2R^3Ry)] 

Poisson's ratio 
atmospheric density 


A 

B 


x,y 


solid A 
solid B 

coordinate system defined in report 



DIMENSIONLESS GROUPING 


The variables resulting from the isothermal EHL point contact analysis (ref. 9) 
are 

(1) Effective radius in the x direction, R^, mm 

(2) Effective radius in the y direction, Ry, mm 

(3) Film thickness, h, mm 

2 

(4) Effective elastic modulus, E', N/mm 

(5) Surface velocity in the x direction, u, mm/sec 

(6) Surface velocity in the y direction, v, mm/ sec 

2 

(7) Atmospheric viscosity, 77^, N* sec/mm 

(8) Viscosity pressure index, Z, dimensionless constant 

(9) Normal applied force, F, N 

(10) , (11) Constants used to define the density of the fluid, a and mm^/N 

It has been found (e. g. , see ref. 10) that density has very little effect on minimum 
film thickness for line contact situations; thus, one may assume the same is true for 
the point contact situation. Even though the compressibility effect is still considered in 
the EHL theory (ref. 9), the constants {a and j3) used to define the fluid in the density 
equation will not be used in the minimum film thickness formulation. Therefore, the 
number of variables shown was reduced to nine. 

In the EHL theory developed in reference 9, the viscosity of the fluid is defined by 
Roelands' (ref. 11) formula. Before Roelands* work the effect of pressure on viscosity 
had been accounted for simply by means of a viscosity-pressure coefficient. Roelands 
points out that the more general solution would be possible by accounting for the 
viscosity-pressure dependence of a given oil by means of its asymptotic pressure 

Piv,as 


Tv, as 



V 


Also, Blok (ref. 12) in 1961 arrived at the very important conclusion that all EHL 
results achieved hitherto for an exponential viscosity-pressure dependence can, to a 
fair approximation, be generalized for any given nonexponential dependence simply by 
substituting the reciprocal of the asymptotic isoviscous pressure Vp^y for th® 
viscosity pressure coefficient occurring in those results. From the nine variables the 
following six dimensionless groupings can be written: 

(1) Dimensionless film thickness: 



Rx 


( 1 ) 


4 


where 


— = ^ 

^Ax ^Bx 

The radii of curvature defined in equation (2) are shown in figure 1. It is assumed that 
convex surfaces (fig. 1) have positive curvatures and that concave surfaces have negative 
curvatures. Figure 2 shows the various components of the film thickness along with the 
central and minimum film thicknesses (H^ and As can be seen from this fig- 
ure, Hp and differ, so we will be interested in defining both. 

(2) Dimensionless load parameter: 


W = 


F 

E'R^ 


where 


E’ = 


1 - 


E, 


1 - IV 


B 


E 


B 


(3) 


(4) 


(3) Dimensionless speed parameter: 


U = (5) 

E*R^ 

where 

V = + v^ (6) 

(4) Dimensionless material parameter: 

G = (7) 

Piv, as 

where pj^^ is the asymptotic isoviscous pressure obtained from Roelands (ref. 11). 
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(5) Spin to roll ratio: 


6 = tan"^^-^ (8) 

(6) Ellipticity parameter: 

k = ^ (9) 

b 

The ellipticity parameter is determined entirely from the definition of the radii of cur- 
vature (r» . rT^ , r* , r„ ), and the derivation can be found in reference 9. Therefore, 
'Ax’ Bx’ Ay’ By” 

the dimensionless film thickness equation can be written as 

H = f(W, U, G, e, k) (10) 

The most important practical aspect of the theory developed in reference 9 is the deter- 
mination of the minimum film thickness in a EHL point contact. Therefore, it is impor- 
tant to observe how the minimum film thickness will vary as these dimensionless group- 
ings are varied. 


NUMERICAL SOLUTIONS 

The objective of the work presented in this report was to keep W, U, G, and 9 
constant in equation (10) while varyii^ k and to observe the effect on the minimum film 
thicknesses. The values of the fixed dimensionless parameters are 


U = 0. 1683x10"^^ 

(11) 

W = 0.3686x10"® 

(12) 

G = 4522 

(13) 

0-0 

(14) 


The variables and their values that were common to all the results are shown in table I. 
In this table the values of t?q, Z, and p.^ correspond to those of a paraffinic min- 
eral oil. Also, the velocity in the y direction is zero; therefore, pure rolling exists. 
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To keep the dimensionless parameters U, W, G, and Q constant while varying 
the ellipticity parameter, it was necessary for the radius of curvature r^y to vary as 
shown in table II. Note that the radius of curvature is negative thereby implying the 
solid to be concave. Also shown in table II is the maximum shear stress. 


RESULTS 

The theoretical analysis presented by the authors in an earlier report (ref. 9) was 
programmed on the digital computer along with the input variables. A contour plotting 
program was also developed. 

Figures 3 and 4 show the contour plots of the dimensionless pressure and film thick- 
ness for ellipticity parameters k of 2, 4, and 6. The + symbol indicates the center of 
the Hertzian circle. Note that, because of dimensionless representation of the X and 
Y coordinates, the actual Hertzian contact ellipse becomes a dimensionless Hertzian 
circle regardless of the ellipticity parameter. The Hertzian contact circle is denoted on 
each figure by asterisks. In the upper left corner of each figure the contour label and 
its corresponding value are given. 

In figure 3(a) for k = 2 a pressure spike emerges. This is denoted by contour 
label c at the back end of the contact. In figure 3(b) for k = 4 the pressure spike oc- 
curs at two locations an equal distance from the axial centerline. In figure 3(c) no pres- 
sure spike is evident for the contour values choosen. 

As the ellipticity parameter is decreased (going from fig. 4(c) to (b) to (a)), the 
minimum film thickness contours move away from the axial center, and two minimum 
areas appear. These results produce all the essential features of the previously re- 
ported experimental observations based upon optical interferometry. In particular, the 
two "side lobes" in which minimum film thickness region occur are shown to emerge in 
the theoretical solutions. 

The effect of ellipticity parameter k on the dimensionless film thickness is shown 
in figure 5. The physical meanings of the minimum film thicknesses H^ and H^^^ 
are shown in figure 2. The dimensionless central film thickness H^ was obtained from 
Archard and Cowking (ref. 7). Writing the Archard-Cowking formula in terms of the 
dimensionless grouping of this report, one obtains 

H^ = 2. 04(AGU)®- (15) 
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where 


A = 


1 


1 + 



(16) 


The Archard-Cowkii^ formula tends to underestimate the central film thickness 
by up to 25 percent (see fig. 5). It also tends to underestimate the minimum film thick- 
ness for values of k greater than three and to greatly overestimate the minimum film 
thickness for k’s less than three. The largest discrepancy occurs for k = 1 where 
the difference between the Archard-Cowking formula and the results of this report is 
over 400 percent. It should be remembered, however, that the Archard-Cowking for- 
mula predicts the central film thickness. 

Also shown in figure 5 is the Dowson-Higginson line contact results where the min- 
imum film thickness is written as 


H 


D 


1 . 6 


‘W 


0. 13 
D 


(17) 


where 



and 


F = 


Force 
Unit length 


If the unit length is taken as the major axis of the contact ellipse. 


F = X 

2a 


(18) 


(19a) 


(19b) 


From figure 5 we see that the minimum film thickness results of the present paper does 
approach the Dowson-Higginson asymptotic value. 
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FILM THICKNESS FORMULAS 


Having determined how the dimensionless fihn thickness varies with the ellipticity 

parameter for a given U, W, G, and 6, the next task was to determine what a more 

general formulation for the minimum film thickness might be. Table HI shows the effect 

of the ellipticity parameter on the new minimum film thickness formula (eq. (21)). In 

this table a film thickness is introduced that has not appeared before: H . t which is 

min, L 

the dimensionless minimum film thickness for line contact obtained from equation (17) 
and which was found to be 6. 955x10"®. 

In table HI a percent deviation of the minimum film thickness is given where 


min 


1.6 e"®-®^ ^ 


^min, L/ / V^min, L/ 


( 20 ) 


The value in table III are observed to be small. Therefore, from this table gen- 

eral formulas for the dimensionless minimum film thickness of point contacts can be 
written as 


“min = “min, L<1 ' 1- « e'"' ■'j 


( 21 ) 


Another approach to determine these formulas would be to use the radius of curva- 
ture ratio Ry/R^ instead of the ellipticity parameter k. The results are tabulated in 

table IV where is defined by 


P . = 
min 




2/3 


H 


min 


H. 


min, Lj 


( 100 ) 



min 


min, L/ 


The value of in table IV is small. It was found that by letting 


( 22 ) 


k = 



(23) 


in equation (21) it did not greatly affect the value of in table IV. Therefore, in 

terms of the radius of curvature ratio R /R , the dimensionless minimum film thick- 

y ^ 

ness can be written as 
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-0. 62(R /R )2/3' 


H ““ H 

min min, L 

1 - 1. 6 e y ^ 

(24) 


In reference 6 Bahadoran and Gohar describe the use of optical interferometry in 
film thickness experiments. In these experiments the radius of curvature ratio 
of the contacting solids was varied. These investigators in their minimum film thick- 
ness measurements used radius of curvature ratios (Ry/R^^) varyir^ from 6. 67 to 22. 66. 
Over this range they make the following conclusion: ”the minimum film thickness 
measured is seen to depend hardly at all on Ry/R^j which charges by a factor of more 
than 3. ” Using equation (24), table V was constructed where the effect of radius of cur- 
vature on percent difference in minimum film thickness for line contact relative to that 
of point contact is shown. From this table it is seen that for Ry/R^ = 6. 67 and 22. 66 
the minimum film thickness varied by less than 17 percent. Therefore, for 
Ry/Rjj 2 : 6. 67 the findings of the present work are in complete agreement with refer- 
ence 6. However, for Ry/R^ = 1 the minimum film thickness differs by 86 percent 
from the line contact solution. One might speculate that, if the investigators of refer- 
ence 6 had used smaller values of radius of curvature ratios (Ry/R^^), they would have 
then seen larger differences in the minimum film thickness measured. The same sort 
of conclusions can be derived from figure 5; that is, for k > 3 (Ry/R^^. > 5. 2) there is 

little change in the minimum film thickness. However, for k < 3 (R /R < 5. 2) there 

y ^ 

is a substantial change in minimum film thickness. 


SUMMARY OF RESULTS 

A numerical solution of the isothermal elastohydrodynamic problem for point con- 
tacts presented herein reproduces all the essential features of the previously reported 
experimental observations based on optical interferometry. In particular, the two side 
lobes in which minimum film thickness regions occur emerge in the theoretical solu- 
tions. 

The influence of the ellipticity parameter on solutions to the point contact problem 
has been explored. The ellipticity parameter k was varied from 1 (a ball on a plate) to 
8 (a configuration approachii^ line contact). It was shown that the minimum film thick- 
nesses can be related to the well known line contact solutions by remarkably simple ex- 
pressions involving either k or the effective radius of curvature ratio Ry/R^- R the 
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II 


additional subscript L is used to denote the line- contact solution, it is found that the 
minimvun film thickness H . for a point contact is 


H = H 
min "min, L 


1 - 1. 6 e 


-0. 62 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, November 17, 1975 
505-04. 
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TABLE I. - VARIABLES COMMON TO ALL RESULTS PRESENTED 


Modulus of elasticity of solid A, E^, N/mm^ 

2X10® 

Modulus of elasticity of solid B, Eg, N/mm^ 

2X10® 

Poisson's ratio of solid A, dimensionless 

0.3 

Poisson's ratio of solid B, i>g, dimensionless 

0.3 

Radius of curvature of solid A in x direction, r^, mm 

11. 113 

Radius of curvature of solid A in y direction, r . , mm 

11. 113 

Radius of curvature of solid B in x direction, rg^, mm 

OO 

Surface velocity in y direction, v, mm/sec 

0 

Lubricant viscosity at atmospheric condition, N- sec/mm^ 

0.411x10*'^ 

Viscosity pressure index, z, dimensionless 

0.67 

Constant used in defining density, <j, mm^/N 

5. 83x10“^ 

Constant used in defining density, /3, mm^/N 

1. 68x10*® 

Asymptotic isoviscous pressure, pj^ N/mm^ 

45. 22 

Normal applied force, F, N 

10 

Surface velocity in x direction, p, mm/sec 

100 


TABLE II. - VALUE OF VARIABLES AS 
ELLIPTICITY PARAMETER IS VARIED 


[Dimensionless parameters: speed, U, 
0. 1683X10'^^; load, W, 0.3686x10*®; 
material, G, 4522.] 


Ellipticity pa- 
rameter, 
k, 

dimensionless 

Radius of curva- 
ture of solid B in 
y direction, 

'’By' 

mm 

Maximum shear 
stress, 

*^max’ 

N/mm® 

1 

—OO 

574. 03 

2 

-17. 143 

417.32 

3 

-13.698 

353. 85 

4 

-12.634 

316. 96 

6 

-11.861 

278.31 

8 

-11. 570 

246. 90 



TABLE m. - EFFECT OF ELLIPTICITY PARAMETER ON 


MINIMUM FILM THICKNESS FORMULA 


[Hmin,L = 6. 955x10-6.] 


Ellipticity 

parameter, 

k 

Minimum film 
thickness 
ratio, 

M /tj 

min' min, L 

1 - 1. 6 e-°- 62 k 
(a) 

Percent de- 
viation of 
minimum 
film thick- 
ness, 

p 

•^m’ 

percent 

1 

0. 1366 

0. 1393 

1. 968 

2 

.5320 - 

. 5370 

.937 

3 

.7692 

. 7509 

-2.376 

4 

.8805 

.8660 

-1. 646 

6 

.9449 

. 9612 

1. 728 

8 

. 9777 

. 9888 

1. 133 


^See eq. (20). 


TABLE IV. - EFFECT OF RADIUS OF CURVATURE RATIO ON 
MINIMUM FILM THICKNESS FORMULA 


T = 6. 955x10'^.] 
^ min, Jj 


Radius of 
curvature 
ratio. 

Minimum film 
thickness 
ratio, 

•TT /tj 

min min, L 

-0. 62(R /R,^)^/^ 
1 - 1. 6 e y ^ 

(a) 

Percent de- 
viation of 
minimum 
film thick- 
ness, 

p 

^m’ 

percent 

1.000 

0. 1366 

0. 1393 

1. 968 

2.843 

. 5320 

. 5390 

1.316 

5.301 

. 7692 

.7571 

-1. 573 

8.303 

.8805 

.8741 

-. 727 

15.85 

.9449 

. 9680 

2. 445 

25. 29 

.9777 

.9923 

1. 493 


^See eq. (22). 
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TABLE V. - EFFECT OF RADIUS OF CURVATURE 


RATIO ON PERCENT DIFFERENCE OF 
MINIMUM FILM THICKNESS FOR 
POINT CONTACT RELATIVE TO 


THAT OF LINE CONTACT 


Radius of cur- 
vature ratio, 
Ry/R, 

M _ 1100 

Elliplicity pa- 
rameter, 
k = (Ry/R^)^/^ 

\ “min, L/ 
percent 

1 

86.07 

1 

6.67 

17. 78 

3. 543 

22. 66 

1. 117 

8. 008 


F 



Figure 1. - Geometry of contacting elastic soiids. 


DEFORMED 



thickness and central and minimum film 
thicknesses. 






CONTOUR LABEL AND VALUE 
A = 0.0000040 
B = 0 0000042 
C =■ 0 0000046 
D = 0. 0000052 
E = 0 0000060 
F = 0 0000074 
G = 0 0000090 





(a) Ellipticity parameter, k. 2. 


CONTOUR LABEL AND VALUE 
A = 0 0000063 
B = 0. 0000066 
C = 0 0000070 
D - 0, 0000075 
E = 0. 0000083 
F = 0. 0000093 
G = 0. 0000110 


CONTOUR LABEL AND VALUE 
A = 0. 0000067 
B = 0. 0000070 
C = 0. 0000073 
D = 0. 0000078 
E = 0. 0000086 
F = 0. 0000100 
G = 0. 0000120 


(b) Ellipticity parameter, k, 4. 


(c) Ellipticity parameter, k, 6. 


Figure 4. - Contour plot of dimensionless film thickness (H = h/Rj,). 



ELLIPTICITY PARAMETER, k = a/b 

Figure 5. - Effect of ellipticity parameter on dimensionless minimum film 
thickness. 
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